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ABSTRACT 

With the advent of optical interferometers that will be coming online in the 
next decade, radial velocity searches for extra-solar planets will be complemented 
by high angular resolution astrometric measurements. In this paper, we explore 
the sensitivity of the astrometric technique, and develop an analytical under- 
standing of the sensitivity in the regime where the orbital period is longer than 
the time-baseline of the survey. As in an earlier paper that dealt with the radial 
velocity technique, our analysis of the astrometric technique utilizes both the 
orbital amplitude and phase. 



1. Introduction 



Radial velocity surveys of nearby stars have produced spectacular discoveries of extra- 
solar planets over the past two decades (see Marcy, Cochran & Mayor 2000). As these 
searches continue into the next decade, they will be supplemented by high angular resolu- 
tion searches for the astrometric signatures of planetary orbits (e.g. Banner & Unwin 1999; 
van Belle et al. 1998; Horner et al. 1999). While radial velocity (RV) technique works by 
detecting the reflex acceleration of the star due to the orbital motion of the companion, 
astrometric technique detects the apparent reflex motion of the star on the sky. In a pre- 
vious paper (Eisner & Kulkarni 2001, hereafter EK2001), we examined the RV technique 
in detail, paying particular attention to the regime where the time-baseline of the obser- 
vations is shorter than the orbital period of the planet. This regime is interesting because 
while the longest baseline of current RV surveys is ~ 12 years, one expects giant planets to 
form in the colder regions of the proto-planetary nebula, and thus we expect such objects 
to possess periods of many years to centuries (Boss 1995). While RV technique is most 
sensitive to short-period, close-in planets, astrometric technique obtains optimal sensitivity 
to long-period massive planets expected from current theories. Moreover, in EK2001, we 
demonstrated that one can achieve a significant improvement in sensitivity (over current 
methods) if the orbital amplitude and phase are included in the analysis. 
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In this paper, we develop an analytical understanding of the sensitivity of astrometric 
technique in this so-called "long-period regime", and address the issues of detection and 
detection efficiency. To the best of our knowledge, the last authors who dealt with the sen- 
sitivity of the astrometric method were Black & Scargle (1982), who used a periodogram 
approach (Black & Scargle 1982). As in EK2001, we will argue below that a least-squares ap- 
proach is preferable to the periodogram. Essentially, the least-squares approach (in contrast 
to the periodogram method) offers the most general approach, and requires no modification 
in the long-period regime. 

The plan of the paper is as follows. First, wc address the case of edge-on orbits in §2. In 
§2.1, we summarize the basic equations and describe the least-squares approach in §2.2. In 
§2.3 we provide analytical estimates for obtaining a detection in the absence of any signal, 
and in §2.4 we obtain estimates for minimum detectable signals in the presence of noise. We 
repeat this analysis for the case of face-on orbits in §3 and for orbits with general inclination 
angles in §4. We conclude in §5. 



2. Edge-On Orbits 

2.1. Basic Equations 

We will assume circular orbits throughout this discussion, and we will neglect the effect 
of annual parallax. The astrometric signature of an edge-on circular orbit is given by 

e{t)=A sin + +Xt + ii, (1) 

where t is the time; A and /j, are the proper motion and position offset of the planetary 
system, respectively; and 

^~ D U^r'^^V ■ ^ ' 

Here, D is the distance to the system, is the mass of the star, and Mp is the mass of 
the planet. Thus, we can express the sensitivity (defined as the minimum-mass planet that 
can be detected) of the astrometric technique in terms of the amplitude of the astrometric 
orbital signature: 

M, = Da[-^) . (3) 

However, it is more difficult to identify planets with long periods than Equation 3 might 
suggest. In the so-called "long-period regime", the total baseline of the observations is less 
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than the orbital period of the planet, and thus we only observe some fraction of the full 
orbital amplitude. Thus the sensitivity is expected to critically depend on the orbital phase. 
As we will show explicitly below, astrometric technique (applied to an edge-on orbit) will be 
the most sensitive when the orbital curvature is maximal (i.e. = [n + l/2]7r). In contrast, 
when the orbital phase is close to rnr (i.e. a pure sine wave), the lack of orbital curvature 
makes it difficult to distinguish the signal from the unknown proper motion of the system. 



The signal analysis for the astrometric technique consists of fitting the observations 
to the model specified in Equation 1. As noted by several authors (e.g. Scargle 1982; 
Nelson & Angel 1998; Eisner & Kulkarni 2001) the most optimal fitting is obtained by using 
the technique of least-squares. First, we derive a linear model equation from Equation 1; 



Here, Ac = ^sin(/). As = ^cos0. Black & Scargle (1982) estimated fi and A first and 
then applied periodogram analysis to the residual data (i.e. with contributions from /x and 
A subtracted). However, as discussed in EK2001, it is necessary to determine all of the 
parameters simultaneously, especially in the long-period regime. 

Using Equation 4 as our physical model, we perform the following analysis to search for 
a signal in a data set. The parameters of the data set are as follows: the duration of the 
survey is Tq, 6'{tj) is the measured astrometric offset at epoch tj, and Uq is the number of 
measured epochs. With no loss of generality, we let the time go from t — — 7o/2 to To/2. To 
find the four unknowns, Ac, As, A, and we minimize 



with respect to Ac-, As-, A, and ji. Here, uj — '^'k/t., where r is the orbital period. This yields 
a matrix equation for the four unknowns, which must be solved numerically in the general 
case. However, in the short-period regime where we observe many orbital cycles, many of 
the summations will average to zero, and this matrix equation will become approximately 
diagonal. In this case, the fit parameters are given by 



2.2. Least Squares Fitting of Sinusoids 



6{t) — AcCOs{ujt) -\- As s\n{ujt) -\- Xt -\- ji. 



(4) 



no — 1 



= ^ [9'{tj) - AcCOs{ijjtj) - AsSm{ijjtj) - Xtj - //]^ 

j=0 



(5) 




(6) 
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As = — ^ 9'{tj) sm{(jjtj), 



A 
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no 



(7) 
(8) 
(9) 



j=0 



2.3. Type I Errors 

Type I errors describe the probability that a high amphtude will be obtained even when 
no signal is present in the data. We assess the statistics of Type I errors, assuming that our 
data set consists of Gaussian noise with a mean of zero and a standard deviation of ctq. 



2.3.1. Short-Period Regime 

Prom Equations 6-9, we see that Ac, As, A, and /j, are sums of Gaussian variables, and 
thus by the Gaussian addition theorem, the four derived parameters must also be Gaussian. 
Specifically, Ac and As obey a Gaussian distribution with a mean of zero and a variance of 

a' = ^, (10) 
no 

where no is the number of measurements. Denoting by Ais the value of \Ac\ (or \As\) that 
is exceeded in 1% of cases, wc note that 

Ais = 2.61(7 = 3.69 no ^^Vq. (11) 

From Equation 9, we see that /i obeys a Gaussian distribution with zero mean and a 
variance of a^/no- Denoting by the value of that is exceeded in 1% of cases, 

fxis = 2.61(7 = 2.61 no ^^Vq. (12) 

Since the statistics of A are somewhat more complicated, we will infer the approximate 
value of Ai5 from the value of /iis'. 

Ai, « ^ = 5.22 aon-'^^To\ (13) 
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2.3.2. Long-Period Regime 

In this regime, the analytical expressions for Ac,As, A, and fi that we derived above 
no longer hold, and so we resort to simulations. To this end we create a synthetic data set 
consisting of Gaussian noise with ctq = 100 /xas (we simulate this measurement accuracy 
because future instruments such as FAME and Keck Interferometer will probably attain 
similar accuracies; Horner et al. 1999; van Belle et al. 1998). We sample the synthetic data 
at one month intervals for Tq — 12 years. We explore orbital periods from 5 years to 100 
years, choosing the interval between sampled periods so as to result in a one radian decrease 
in the total number of orbital cycles over the length of the observations, 

For each period, r, we simulate = 1000 data sets and carry out the least-squares fit. We 
set Aci equal to the 10**^ highest value that arises. Clearly, 99% of fitted Ac will lie below 
this value. Using the same method, we also determine the 99% confidence level for the other 
fitted parameters, which we denote by Asi, Ai, and /xi. 

We find that the values of the fitted parameters are in excellent agreement with the val- 
ues predicted by Equations 6-9 in the short-period regime (Figure 1). In order to understand 
the behaviors of the fitted parameters in the long-period regime, we begin by examining how 
the covariance of the fitted parameters depends on period. We expect that in the short- 
period regime, all four parameters will be completely uncorrelated, since we observe many 
orbital periods and can thus easily distinguish the sinusoidal motion from the linear proper 
motion term and the constant offset term. This expectation is verified by the simulations 
(Figure 2). 

In the long-period regime, in contrast, strong correlations may arise between the fitted 
parameters. From Figure 2, we conclude the following for the long-period regime: 

1. Ac and As are correlated, and select orbital phases close to 0°. 

2. Ac and /i are anticorrelated. 

3. As and A are anticorrelated. 

4. Ac and A are uncorrelated (not shown in figure). 

5. As and /i are uncorrelated (not shown in figure). 

Previous authors noted that these correlations would arise, and that the orbital phase would 
become important in the long-period regime (Black & Scargle 1982). However, neither they 
nor others have explored the full- implications of this fact. Below, we look into this issue in 
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more detail, and develop an analytical understanding of how the fitted parameters depend 
upon orbital period. 

We now provide a physical explanation of the results displayed in Figure 2. When 
T <^ Tq, sinusoids with random phase can be fit to Gaussian data, where the amplitude of 
the sinusoids is constrained by the vertical scale of the data (~ ctq). However, when r > Tq, 
this is no longer necessarily the case. The maximum amplitude is obtained when a sine wave 
centered around t = is fitted. This is because for r > To, sin(7rTo/T) ~ ttTq/t, which can 
be completely absorbed into the linear proper motion term. 

Specifically, the maximum amplitude of the sine wave that can be fitted is limited by 
how much sin(27ri/T) deviates from a straight hue in the range from t — [0, ro/2]. This 
deviation is given by 

_ sm — . (15) 



We note that the maximal deviation occurs when r > 3To/4, and thus the long-period regime 
"begins" where r = STq/A. The largest possible amphtude of the fitted sine wave (As) in 
the long-period regime is that for which this deviation is approximately equal to the vertical 
scale of the noise (~ ctq). Since Equation 15 tells us the fractional observed amplitude, the 
real amplitude of the fitted sine is 

ifr„/T-sm(iTT„/T) ^ ' 

Here, Agi is the value of l^^l that is exceeded in 1% of least-squares fits to Gaussian noise, 
Ais is the corresponding value of Agi in the short-period regime (Equation 11), and the 
specific normalization is chosen so that at the boundary of the short and long-period regimes 
{t — 3To/4), Agi is continuous. 

We can understand the behavior of Ac in the long-period regime using a similar analysis. 
Specifically, the maximum amplitude of a cosine wave centered around t = is constrained 
by how much this cosine deviates from a constant in the interval [0,To/2]. Following the 
same lines as the analysis for As, we find that this amplitude is given by 

^1 = 1 — ^4Vtt- (1^) 

1 - COS (ttTo/t) 

Here, Ad is the value of \Ac\ that is exceeded in 1% of least-squares fits to Gaussian noise, 
Ais is the corresponding value of Ad in the short-period regime (Equation 11), and 2 is a nor- 
malization factor (chosen to enforce continuity between the short and long-period regimes). 
We note that in the long-period regime sin(a;) — > x faster than cos(a;) — >• 1, and thus the 
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covariance of the sinusoidal signal with the constant offset term is less significant than the 
covariance with the proper motion term. 

In light of these analytical expressions for Ac and As, it is not difficult to find expressions 
for A and /i. In particular, because the input data consists of Gaussian noise with a mean 
of zero, we know that A and /i must be chosen so that the resultant fit is also, in the mean, 
zero. The anti-correlation of Ac and /i and that of As and A confirm this expectation. 

We must choose a A that "cancels" the fitted sine- wave component. In the long-period 
regime, this condition can be expressed as AsSm{u!t) ~ AsOut ~ Xt, or A ~ As/t. How- 
ever, because the data has some vertical scale, we must vertically "re-center" the data after 
subtracting As/t. Since the total vertical scale is ~ 2(To, we must add (Jq after subtracting 
As/t. Thus, we expect that for r > 3To/4, 

"4t \ irr„/T-sin(xTo/T) J ^ ' 

Similarly, is chosen to "cancel" the cosine component. Thus, ^ is the negative co- 
sine component (given by Ac) plus (Tq (where ctq is needed to reproduce the data, which is 
vertically centered around 0). Thus, for r > Tq, 



Ail = A*is 



1 - cos (ttTo/t) 



- 1 



(19) 



Equipped with the behaviors of A,.i and Agi, we now know the analytical behavior of 
the Type I errors in the long-period regime. Specifically, the error ellipse within which 99% 
of fitted signals lie, ei, is completely described: Agi gives the semi-major axis, and Ad gives 
the semi-minor axis. 

Amplitude— Phase Analysis. The standard signal analysis usually focuses on the statistics 
of the orbital amplitude A = {Al+AlY^"^. This is adequate when Ac and As arc uncorrelatcd 
(e.g. in the short-period regime with adequate sampling). However, as can be seen in Figure 
2, the Ac- As distribution becomes highly elliptical in the long-period regime. Thus in this 
regime it pays to undertake an analysis that utilizes amplitude and phase. The fact that 
the amplitude-phase ellipse (Figure 2) is significantly smaller than the amplitude-only circle 
in the long-period regime clearly indicates the superiority of an amplitude-phase approach. 
Therefore, we will use an amplitude-phase approach in the remainder of this discussion. See 
EK2001 for a more detailed comparison of amplitude-only and amplitude-phase approaches. 
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2.4. Type II Errors 

In §2.3, we computed the probability of detecting an apparent signal generated purely 
by noise (i.e. Type I probabilities). We now consider Type II probabilities- the probability 
of failing to detect a genuine signal due to contamination by noise. To this end, we simulate 
a data set that consists of signal and noise: 

e'{ti) = Asm{ujti + N{ti), (20) 

where A is the amplitude of the signal, (f) is the phase, and N{ti) is the Gaussian noise. We 
let be drawn from a uniform distribution in the interval [0, 2ti] , an appropriate assumption 
for circular orbits. We choose an initial signal amphtude of cro/2, and then do = 1000 
least-squares fits (with the same parameters as in §2.3). 

In §2.3, we showed that in the long-period regime, the Ac-As distribution is elliptical. 
In view of this, we define the Type I errors by an ellipse ei; this ellipse contains 99% of 
the zero-signal simulations (§2.3.2). In our analysis of Type II probabihties, we consider all 
least-squares fits (of the signal specified in Equation 20 to the model specified by Equation 
4) that lie within ei to be indistinguishable from those produced by noise. We increment the 
input signal amplitude until a certain fraction of fitted points lie outside of ei. For example, 
the signal amplitude at which 99% of fitted points lie outside of ei is denoted by ^499. Plots 
of Agg, Ago, and A50 (as a function of orbital period) are shown in Figure 3. 



3. Face-On Orbits 

In the case of face-on orbits, we have twice as many measurements as in the edge-on 
case, since there are two orbital dimensions. Thus, our model will contain both the x and y 
motion of the orbit: 

9xit) = Ac cos{ut) + As sin(cjt) + \xt + i^x, (21) 

6y{t) = -As COS{U}t) + Ac Sm{u}t) + Xyt + jly. (22) 
— * 

As above, Ac — Asmcf), As — Acoscf), A = [A^;, Aj^] is the proper motion, and /2 = [//j,, ji^ is 
the position offset. 

Using the same parameters as in §2.2, we determine the unknown model parameters by 
minimizing 

no-l 

X"^ = '^{[O'xitj) - Ac cos{utj) + As sm{utj) - Xxtj - 

j=0 
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[Oy{tj) + As cos(u;tj) — Ac sm(u;tj) — Xytj — //J^} (23) 
with respect to Ac, As, X^, \y, i^x, and /ly. 

3.1. Type I Errors 

Because a face-on orbit has two observable dimensions, the analysis of Type I errors will 
be somewhat more comphcated than for the face-on case. First, we note that in the short- 
period regime where r < Tq, the least squares fit essentially has twice as many measurements 
to work with. Thus, we expect that Ais will be 2~^/^ times its value for the edge-on case: 

Ais = 2.61(7 = 2-^/^ X 3.69 tIq ^^Vq. (24) 

However, since A and fi must be determined for each orbital dimension independently, /^i^ 
and Ais will still be given by Equations 12 and 13. 

When T > Tq, the full amplitude of a cosine centered around t = cannot be observed. 
However, in the case of face-on orbits, we can observe the full amplitude of the orbit in the 
orthogonal direction, since it is phase-shifted by 90°. The full amplitude of this sine wave 
can be observed until r > 4To/3 (to see this, find the value of r for which cos(7rTo/r) = 
-sin(7rTo/r)). 

For face-on orbits, the long-period regime begins when r > 4To/3, since it is here that 
the full amplitude of a sinusoid cannot be observed in cither orbital dimension. Moreover, 
since the two dimensions are phase-shifted by 90°, we expect the expressions for Ad and Agi 
to be the same. Since sin(a;) x faster than cos(x) — > 1, the deviation of a cosine wave 
from unity will be the important constraint on the possible fitted amplitude. Thus, Ad and 
Asi are both given by Equation 17 in the long-period regime. 

A potential complication is that at some point, we are gaining virtually no new infor- 
mation from one of the orbital dimensions (since the constraint on the size of a cosine is 
so much stronger than that on the size of a sine wave). So, we expect that at some pe- 
riod, our sensitivity will drop roughly by a factor of 2^/^, to reflect the fact that half of the 
measurements become useless. We model this behavior as follows: 

{Ais for r < To 

V2Au for To < r < 4To/3 (25) 
1-^^/.) forr>4To/3 
As illustrated by Figure 4, this function is a good fit to the simulations. 

Since Asi now depends on the maximum amplitude of a cosine wave, rather than the 
amplitude of a sine wave, we must modify our expression for Ai accordingly. Specifically, 
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while the slope of a sine wave centered around t — goes approximately as cut, the slope of 
a cosine wave goes as cu^t^. Thus, we expect that the expression for Ai in the case of face-on 
orbits will differ from that for edge-on orbits by a factor of cu. Specifically, we find that for 

T>To, 

A:^aJ^)72 (26) 

\T J [ ttTq/t - sm (vrTo/r) J 

This analytic function is shown along with the simulated data in Figure 4. 

We expect that fii will be the exactly the same as in the case of edge-on orbits, since 
the behavior of Ad is unchanged. This expectation is confirmed by Figure 4. 

From these analytic expressions for Ac and As, we see that the sensitivity of astrometric 
technique to face-on orbits is significantly better than that for edge-on orbits. This is well 
illustrated through a comparison of Figures 2 and 5. As discussed above, the Type 1 error 
ellipse, ei, is defined as the ellipse in the Ac-Ag plane for which 99% of fitted parameters 
lie inside. In the case of edge-on orbits, this ellipse is very eccentric, since the semi- major 
axis is given by Equation 16 and the semi-minor axis is given by Equation 17 (Figure 2). In 
contrast, in the case of face-on orbits, ei is a circle whose radius is given by Equation 25. 
Thus, for face-on orbits, we are equally sensitive to signals of arbitrary orbital phase. 



3.2. Type II Errors 

As in §2.4, we compute the probability that a real signal will not be detected due to 
contamination by noise. For a face-on orbit, we simulate the following data set: 

e'^iti) = Asm{u;ti + + N^iU), (27) 

e'yiU) = AcosiuU + 0) + NyiU), (28) 

Here, A is the amplitude of the signal, is the (randomly-distributed) orbital phase, and 
N{ti) is the Gaussian noise. Since the sensitivity of astrometric technique to face-on orbits is 
better than for edge-on orbits by a factor of 2^/^ (in the short-period regime), we choose an 
initial signal amplitude of (To/2^/^. We then do N — 1000 least-squares fits, with the same 
parameters as in §3.1. 

Following our analysis of Type II errors in §2.4, we consider all least-squares fits (of the 
signal specified in Equations 27 and 28 to the model given by Equations 21 and 22) that 
lie within ei to be indistinguishable from noise. The results of this analysis show that the 
sensitivity of the astrometric technique is much better for face-on orbits than for edge-on 
orbits (Figure 6). 
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4. Inclined Orbits 

In the case of generally inclined orbits, we observe the full amplitude in one orbital 
dimension, but only sini x ^ in the other dimension (where i is the orbital inclination 
angle). Labeling the dimension in which we measure the full amplitude by x, our model 
becomes 

9x{t) — AcCOs{u;t) + AsSm(ujt) + Xxt + jix, (29) 

9y{t) = — sin i As cos{ujt) + sin i Ac sin(cjt) + Xyt + fiy. (30) 

As above, Ac = Asm4>, As = Acos4>, X is the proper motion, and jl is the position offset. 

Using the same parameters for our data set as in §2.2, we determine the unknown model 
parameters by minimizing 

no — 1 

X'^ = '^{[O'xitj) - AcCOs{ujtj) - AsSm{ujtj) - Xxtj - 
i=o 

[0'y{tj) + sin i As cos{ujtj) — sin i Ac sm{ujtj) — Xytj — fJ>y\^} (31) 
with respect to Ac, As, Xx, Xy, fix, and fiy. 

4.1. Type I Errors 

In the short-period regime where r < Tq, the least squares fit essentially has twice as 
many measurements with which to work, although half of the measurements are less sensitive. 
Thus, we expect that Type I errors will be somewhere between those for the edge-on and 
face-on cases. However, for simphcity, we will simply use the expressions from the face-on 
case. 

As above. Ad is given by Equation 17 in the long-period regime, since we can always 
observe the full amplitude of a cosine wave in one orbital dimension. However, Asi will be 
different than in the face-on case, since the cosine wave cannot be measured as accurately 
in the other orbital dimension. Thus, we expect that in the long-period regime, Asi will be 
stretched relative to Ad- For orbits with inchnation angles close to zero, there is essentially 
no information from the second orbital dimension, and thus Asi is given by Equation 16. 
However, for inclination angles significantly greater than zero, the information from the 
second orbital dimension becomes important, and Asi is limited by the lower- amplitude 
cosine wave in this second dimension. Specifically, Asi is Ad stretched by 1/ sini. Thus, we 
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have an expression for the aspect ratio of ei as a function of inchnation angle: 



^1 



mm 



sm^ 



(32) 



edge . 



As an example, we determine the Type 1 errors for a 45°-inclined orbit. In this case, 
= Aci/ sini = 2^/^ylci (Figure 7). Since is increased by 2^/"^ from the face-on case, 
Ai must be similarly increased: 



Ai = ^/2Xu 



r 



TT 



smvr 



ttTo/t - sin(7rTo/T) 



(33) 



We expect that /xi will be the exactly the same as in the case of edge-on orbits, since the 
behavior of A^i is unchanged. Analytic functions for the parameters are shown along with 
the simulated data in Figure 8. 



4.2. Type II Errors 

As in §2.4 we compute the probability that a real signal will not be detected due to 
contamination by noise. For inclined orbits, we simulate the following data set: 

e'^iU) = Asm{ujti + 4)) + N^{ti), (34) 

9y{ti) = smiAcos{ujti + 0) + Ny{ti), (35) 

Here, A is the amplitude of the signal, is the (randomly-distributed) orbital phase, and 
N{ti) is the Gaussian noise. We choose an initial signal amplitude, and do = 1000 
least-squares fits, with the same parameters as in §3.2. 

Following our analysis of Type II errors in §2.4, we consider all least-squares fits (of the 
signal specified in Equations 34 and 35 to the model given by Equations 29 and 30) that he 
within €i to be indistinguishable from noise. The results of this analysis for a 45°-inchned 
orbit are shown in Figure 9. 



5. Caveats and Conclusions 

As in our earlier analysis of the sensitivity of the radial velocity technique (EK2001), 
the analysis presented here is based on frequentist statistics. However, the importance of 
the Bayesian approach is gradually being recognized in astronomy. Berger & Delampady 
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(1987) point out that the probabihty of extreme events predicted by the frequentist approach 
(which is what we have essentially computed) is rosier than is warranted. Thus, given a 
mix of stars (with and without planets), the frequentist tail criterion of say 1% (e.g. ei) 
is likely to underestimate the fraction of false detections that arc expected to occur in 
population analyses. These criticisms notwithstanding, we have pursued the frequentist 
approach because of (1) its conceptual simplicity, (2) its ease of use, (3) our own familiarity, 
and (4) our lack of familiarity with the Bayesian approach. 

Wc also acknowledge that out treatment of the statistics was not always completely 
rigorous. Specifically, when one orbital dimension becomes much more sensitive than the 
other, we simply ignore the data from the other dimension. Inclusion of this less sensitive 
dimension in the analysis is likely to lead to a small improvement in the overall sensitivity. 
We feel that this statistical analysis is beyond the scope of this paper, and will not lead to 
any meaningful modifications of the results derived above. 

A final caveat is that in this analysis, we have assumed that the orbital inclination 
angle is known a priori. In general, this is not the case, and the inclination angle must 
be determined simultaneously with the rest of the orbital parameters (using a non-linear 
least-squares fitting method) . While one expects that all of the sensitivities in this paper are 
slightly under-estimated as a result, it is clear that the qualitative behavior of the sensitivity 
in the long-period regime is unaffected. 

Despite the fact that our analysis can be improved upon in some areas, this work provides 
an accurate analytic description of the sensitivity of the astrometric technique in detecting 
outer planets. Moreover, we have understood, both qualitatively and quantitatively, how 
the sensitivity improves as the orbital inclination angle is increased. Of specific importance 
is the fact that in the long-period regime the sensitivity of astrometric technique to face-on 
orbits is more that 2^/^ times better the sensitivity to edge-on orbits, due to the orthogonal 
information supplied by the two orbital dimensions. 

REFERENCES 

Berger, J.O., & Delampady M., 1987, Statistical Sci.,2(3),317 
Black, D.C., & Scargle, J.D., 1982, Ap. J., 263, 854 
Boss, A.P., 1995, Science, 267, 360 

Banner, R., & Unwin, S., 1999, Space Interferometry Mission Book 
Eisner, J.A., & Kulkarni, S.R., 2001, Ap. J., 550, 871 



-14- 



Horner et al., 1999, in lAU Colloquium 176, The Impact of Large-Scale Surveys on Pulsating 
Star Research, eds. L. Szabados & D. Kurtz (ASP Conference Series) in press. 

Marcy, CD., Cochran, W.D., & Mayor, M., 2000, Protostars and Planets IV, eds Mannings, 
v.. Boss, A. P., & Russell, S. S., University of Arizona Press, p. 1285 

Nelson, A.F., & Angel, J.R.P., 1998, Ap. J., 500, 940 

Scargle, J.D., 1982, Ap. J., 263, 835 

van Belle, G.T., Boden, A.F., Colavita, M.M., et al., 1998, SPIE, 3350, 362 



This preprint was prepared with the AAS IM^jX macros v5.0. 



- 15 - 




0.0 0.5 

log(T/To) 



1 .0 





3 


- 




2 


5 - 




2 


- 




1 


5 : 


<^ 

cn 
O 




- 







5 : 







- 




-0 


5 - 





1 


/ ' 




f 






























/ 






/ 











-0.5 0.0 0.5 

log(VTo) 



1 .0 



o 




0.0 0.5 

log(T/To) 



0.0 0.5 

log(T/To) 



1 .0 



Fig. 1. — Illustration of how the fitted parameters depend on orbital period in the case of 
edge-on orbits. The solid lines indicate the 99**^ percentile values of the orbital parameters 
Act Asi A, and yU. Here, we simulated = 1000 data sets with Gaussian noise of zero mean 
and (To = 100 /xas. Data were sampled at one-month intervals for Tq = 12 years. The dashed 
lines show the analytic estimates predicted by Equations 16-19. 
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Fig. 2. — Correlations between fitted parameters for edge-on orbits in the short- and long- 
period regimes. Each dot represents one out of the = 1000 simulations. The duration 
of the astrometric monitoring is Tq — 12 yr. In the short-period regime, r <^ To, all of the 
orbital parameters are uncorrelated. In the long-period regime, strong correlations are seen. 
In particular, orbital phases of approximately 0° are preferred (i.e. the Type I error ellipse 
has a large aspect ratio), Ac and /i are anticorrelated, and As and A are anti-correlated. 
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Fig. 3. — Plots of 7499, ^490, and ^450 and a function of orbital period, for edge-on orbits. ^499 
is the signal amplitude necessary such that least-squares fits to data containing a genuine 
signal plus Gaussian noise will yield detections 99% of the time. We simulated = 1000 
data sets sampled at one-month intervals for Tq = 12 years. The data sets consist of Gaussian 
noise (with zero mean and (Tq = 100 //as) plus signal; the phase of the signal was assumed 
to be randomly and evenly distributed over the range [0, 2tt]. 
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Fig. 4. — Illustration of how the fitted parameters depend on orbital period in the case of 
face-on orbits. The solid lines indicate the 99*'^ percentile values of the orbital parameters 
Ac-, As-, A, and /x, and the dashed lines show the analytic estimates predicted in §3.1. See 
Figure 1 for details of the simulations. 




Fig. 5. — Correlations between fitted parameters for face-on orbits. Each dot represents one 
out of the N — 1000 simulations. The duration of the astrometric monitoring is Tq — 12 yr. 
In the case of face-on orbits, we see that Ac and As remain uncorrelated in the long-period 
regime, and thus the aspect ratio of the Type I error ellipse ei is unity. 
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Fig. 6. — Plots of Agg, Aqq, and A^q as a function of orbital period, in the case of face-on 
orbits. Agg is the signal amplitude necessary such that least-squares fits to data containing 
a genuine signal of random phase plus Gaussian noise will yield detections 99% of the time. 
See Figure 3 for details of the simulations. 




Fig. 7. — Correlations between fitted parameters for 45°-inclined orbits. Each dot represents 
one out of the A'^ = 1000 simulations. The duration of the astrometric monitoring is Tq = 12 
yr. In the case of 45°-inclined orbits, we see that the aspect ratio of the Type I error ellipse 
€i is 2^/^ in the long-period regime. 
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Fig. 8. — Illustration of how the fitted parameters depend on orbital period in the case of 45°- 
inclined orbits. The solid lines indicate the 99*'^ percentile values of the orbital parameters 
Ac-, As-, A, and /x, and the dashed lines show the analytic estimates predicted in §4.1. See 
Figure 1 for details of the simulations. 
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Fig. 9. — Plots of Agg, Ago, and A^q as a function of orbital period, in the case of 45°-inclined 
orbits. Agg is the signal amplitude necessary such that least-squares fits to data containing 
a genuine signal of random phase plus Gaussian noise will yield detections 99% of the time. 
See Figure 3 for details of the simulations. 



